Rocking ratchets in 2D Josephson networks: collective effects and current reversal 
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A detailed numerical study on the directed motion of ac-driven vortices and antivortices in 2D 
Josephson junction arrays (JJA) with an asymmetric periodic pinning potential is reported. Dc- 
voltage rectification shows a strong dependence on vortex density as well as an inversion of the vortex 
flow direction with ac amplitude for a wide range of vortex density around /=l/2 ( f—Ha 2 /$o), 
in good agreement with recent experiments by Shalom and Pastoriza [Phys. Rev. Lett. 94, 177001 
(2005)]. The study of vortex structures, spatial and temporal correlations, and vortex-antivortex 
pairs formation gives insight into a purely collective mechanism behind the current reversal effect. 

PACS numbers: 74.81. Fa, 74.25. Qt, 85.25.Cp 
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Ratchet systems constitute a vast subject of present 
interest, extensively studied from biology to a wide spec- 
trum of fields in physics 0, 0, H, 0, [a]. Very recently, 
promising applications in micro-nano technology have 
motivated a huge amount of additional works, mainly re- 
garding the controlled motion of nanoparticles, colloids, 
electrons or magnetic flux quanta. The main characteris- 
tic of these systems is their ability to transport particles 
with a non-zero drift velocity in the absence of any net 
macroscopic force. The simplest realization of this idea 
is the rectification of an ac-driven single particle motion 
in an asymmetric periodic potential or ratchet potential, 
known as a rocking ratchet. Ratchet effect in supercon- 
ductors (SC) has been, in particular, an attractive sub- 
ject of study since the proposal of using it to remove un- 
desirable trapped vortices producing noise in SC devices 
Q. Notably, several designs of asymmetric pinning po- 
tentials for the motion of flux quanta have been recently 
realized, such as dots-antidots arrays of pinning sites in 
SC films [3]. Even more, ratchet designs without spatial 
asymmetry have been successfully achieved in layered SC 
[§|. Small Josephson junction (J J) systems were inves- 
tigated as well: SQUIDs @], ID ladders 0, [ll|, long JJ 
[pi ] and quasi-lD quantum JJA (l3| . In particular, on 
classical JJA, the first desi gn o f a device for kinks motion 
rectification was proposed [lfj and measured ll| in ID. 

Large 2D JJA, which are excellent systems to study 
statistical mechanics experimentally [14j and where col- 
lective effects could play a relevant role, had not been 
studied until very recently. Rectification experiments on 
2D JJA were reported by Shalom and Pastoriza [HI]. 
They modulated the distance between the SC islands 
thus generating a ratchet potential for vortex motion. 
Dc-voltage rectification was then clearly observed and 
its intriguing features, such as a non trivial dependence 
on both, the applied ac current and vortex density, and 
an inversion of the vortex flow direction, prompt further 
investigation. The aim of this work is to interpret such 
experiments and understand the role of collective effects. 
We show numerically that the overdamped RSJ model 
for large 2D JJA with asymmetrically modulated criti- 
cal currents, reproduces successfully the experiments (l5| 



and allows for a detailed characterization of the device 
useful for performance enhancement and future experi- 
ments. In addition, a complete study of the vortex dy- 
namics gives insight into a non trivial collective mecha- 
nism behind the current reversal (CR) effect. 

We study the overdamped vortex dynamics in square 
2D asymmetrically modulated JJA using the resistively 
shunted junction (RSJ) model [lfj and solving numeri- 
cally the superconducting phase dynamics. The current 
flowing in the junction between two superconducting is- 
lands is the sum of the Josephson supercurrent and the 
normal current: I M (n) = I sin 6> M (n) + 2n % w d %[ n) + 
77^ (n, t), where Iq is the critical current of the junction be- 
tween the two superconducting islands sited in n and n + 
/jona square lattice [n=(n x ,n 2/ ), /x=x,y], M (n)=0(n + 
/i) — 9(n) — ^(n) the gauge invariant phase difference 

with A M (n)=^ J,," ° A-dl, <&o the flux quantum, a the 
lattice parameter, and Rn the normal state resistance. 
Under a magnetic field H, we have A M x yl M (n)= 2tt/, 
/=7Ja 2 /$o being the frustration parameter, a measure 
of the vortex density. The Langevin noise term r/ M models 
the contact with a thermal bath at temperature T and 
satisfies (^(n,t) V (<tO}=^<W<W<Kt-t'). The 
ratchet potential for vortices is introduced by modulating 
the critical currents Io=EjQq/2tt, Ej being the Joseph- 
son coupling energy between islands. We consider a saw- 
tooth modulation Io(n x ) increasing linearly from Ia min 
to Io mam within each period, p, as an approximation to 
the experimental modulation. The ratchet potential am- 
plitude is then measured by AU oc A7o=(/o max — ^o min )- 
The system is driven by I ex t =I ac sin(27ra; ac i) in the y 
direction. Periodic boundary conditions in both direc- 
tions are taken (see model details in fl7j). The dynam- 
ical equations for the superconducting phases are ob- 
tained after considering conservation of the current in 
each node. The resulting set of Langevin equations are 
solved using a 2 nd order Runge-Kutta algorithm with 
time step At=0.1rj (tj=2itcRnIo/®o, where Iq is the 
averaged critical current) and integration time 2.10 6 rj 
after the equilibration transient. We calculate: mean 
dc voltage as (Vdc)~-(dOy(n)/dt) [13] vs both I ac and 
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FIG. 1: Rectified vortex motion: I ac — Vdc characteristics ex- 
amples increasing vortex density, from one vortex in the sam- 
ple (/=1/L 2 ) towards f=l/2. Inset: {V dc ) vs I ac for / '=1/32, 
in log-log scale (right y-axis) and vortex-antivortex pairs den- 
sity normalized by frustration, n va / f (left y-axis). 



/ where <..> stands for both time and space average; 
vorticity at the plaquette n=(n x + 1/2, n y + 1/2), as 
6(n)=— A p x nint[0 Al (n)/27r] [18| (nint[] is the nearest 
integer function). Currents are normalized by Iq, volt- 
ages by RnIo> temperature by Ej/Icb and frequencies by 
(tj) _1 . Most of the results shown are (except indication) 
for p=8, I Om „ i =0.57o, I(w=1.5io, cj qc =10- 4 , T=0.05 
and 32x32 arrays. Other parameters as p=7, 10, 15 
and A/o=0.5, 1.5, 2 were as well investigated. The full 
range of magnetic field was explored, from a single vor- 
tex to the fully frustrated case (f—1/2) in system with 
L=L x =L y =32, 64, 100, 128. 

Results. In FigfT] we show mean dc- voltage, (Vd c ), vs 
the amplitude of the ac current, I ac , for different vor- 
tex densities, /. These results imply rectification: a net 
directional vortex motion parallel to the ratchet modula- 
tion in response to an alternating force. For low /, (Vdc) 
is maximum at an optimal value of I ac and decreases 
slowly for larger I ac as the ratchet potential is effectively 
washed out by the fast vortex motion. This behavior is 
qualitatively the one expected for a single particle rock- 
ing ratchet system Q • Let us note however that even for 
highly diluted vortex systems the intensity of rectification 
is not additive: for example with /=l/32 in a 32 x 32 
array we have 32 vortices but as we can see in FigJTJthe re- 
sponse is ^12 times, and not 32 times the response of one 
single vortex (/=1/L 2 ). This non-additivity is even more 
pronounced for higher /. By increasing frustration from 
/=l/32 to /=l/2 the rectification peak presents both a 
reduction and a shift of the optimal I ac value. A more 
intriguing feature appears when / is further increased. 
For 0.17 < / < 1/2 an inversion of the vortex flow di- 
rection is observed at small I ac values and the window 
with (Vdc) < becomes wider as / is increased. The non 
monotonic behavior observed as a function of / is a clear 
evidence that vortex-vortex interactions are relevant for 
the rectification effect at all /. On the other hand it 
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is well known that the frustrated XY-model, eq 
to our model for 2D JJA under an external H 
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FIG. 2: (a) {V dc ) vs I ac for /=l/4 and /=3/4. (b) Numerical 
{Vdc) contour plots in I ac — f plane at T=0.05Ej /ks and 
cj ac =10~ 4 (l/rj). Red dotted lines represent (Vdc) < — Vo and 
green solid lines, (Vdc) > Vb=5.10~ RnIo- (c) Experimental 
contours at T=3.8K, uj ac —lkHz and Vo=15nV (see [H[). 



also present a finite density of thermally excited vortex- 
antivortex pairs, n va =((\b(n)\) — /)/2, over the magnetic 
field induced vortex pattern. It is then important to 
understand its contribution to the ratchet effect. In the 
inset of Fig.l we plot both (Vdc) and n va /f vs I ac for 
/=l/32. We see that n va increases monotonically and 
saturates to 10% of / at I ac ~ 0.8. We find the same be- 
havior for all / but the saturation value of n va decreases 
with increasing /, being less than 0.01/ for / > 1/8. We 
conclude that the effect of pairs is almost negligible in the 
CR regime, but it is relevant for diluted vortex lattices 
at large applied currents and they must be considered to 
reproduce properly the experimental rectification tails. 

A summary of our numerical results on voltage rec- 
tification vs vortex density is presented in Fig[2] in 
comparison with experiments 15]. In Figl2ja) we 
show numerical curves for /=l/4 and f—1— /=3/4. 
From the symmetries of the frustrated XY Hamilto- 
nian: H(f)—H(—f)—H(l — /) [Hj], we expect symmet- 
rical transport properties around /=l/2. Since vortex 
and antivortex motion is rectified in the same direction 
we have (Vd c )(f)=— (Vd c )(l— /)■ This means that in the 
< / < 1/2 range of magnetic field the dissipation 
arises from the motion of "positive" vortices while the 
flux of "negative" vortices (antivortices) yields dissipa- 
tion in the 1/2 < / < 1 range. In Fig[2^b) contour plots 
of (Vdc) show the occurrence of voltage rectification in the 
Iac — f plane. Current reversals occur only in the regime 
of strong vortex repulsion around /=l/2, in a single in- 
terval, 0.17 < / < 0.5. Comparing Figl^b-c) we see a 
good qualitative agreement with experiments. 

For applications it is also important to accurately 
determine the range of parameters where rectification 
would be observable. In FigEJa) we see a well defined 
saturation plateau for (Vdc) at low frequencies. At high 
Lj ac the fast changes in the potential slope do not al- 
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FIG. 3: Rectification dependences for /=l/8: (a)Onu„ c : at 
7 ac =0.25(*) and J ac =0.75(o) ((V dc ) max , see Fig.l). (b)On 
pinning potential: (Vdc)-I ac vs A/ - Inset: (Vdc)max vs A7 . 



low vortices to adapt to the ratchet potential during a 
cycle. For all I ac the mean voltage thus decays, be- 
coming negligible at w ac >0.0Tg(T/) _1 ~l/T re i. Then we 
predict an observable rectification up to ^lOOkHz in 
good agreement with [IH, if we use typical parameters, 
i?Ar=0.01i7 and Iq=1[iA, to estimate our unit of fre- 
quency (tj) -1 ~30MHz. Rectification can also depend 
on the parameters of the ratchet potential. In FigfS^b) 
we see that (Vd c ) monotonically increase with Alg for 
all I ac , and the optimal rectification increase approxi- 
matively linearly for 1.5 < Alg <2 (see inset) and for 
AIq ~ 2 inversion disappears |19j. We thus see that 
the ratchet effect could be experimentally improved up 
to one order of magnitude by modulating accordingly the 
Josephson couplings. On the other hand we find that rec- 
tification has no appreciable changes for the experimen- 
tal values p=7,10,15 nor sharp commensurability effects 
between / and p. We also analyze the range of temper- 
ature where rectification occurs and size effects. We find 
a wide range of applicability, < T < 0.4, what makes 
JJA advantageous and that the main characteristics of 
our results are almost size independent for L > 32 (lij . 

In order to understand collective transport effects and 
in particular current reversal, we study vortex structures 
and spatial and temporal correlations. In FigfJJa-b) we 
show vorticity snapshots for /= 1/8, 3/4 at I ac = il- 
lustrating the competition between the ratchet poten- 
tial, which induces a modulated vortex density, and 
the repulsive vortex-vortex interaction that favors a 
uniform vortex density. For all / we find that vor- 
tices organize in dense vortex walls strongly pinned at 
the sawtooth minima. Interestingly, when CR starts, 
the static vortex structure present an averaged num- 
ber of vortices in the walls similar to a checkerboard 
pattern, i.e. (n v / waU )^L y /2. The degree of peri- 
odic order is reflected by the intensity of the peaks at 
k=(n27r/p, 0) of the averaged vortex structure factor 
S(k)=(||-^£ fi 6(n)exp(ik-n)| 2 ^ shown in Figgfc). A 
ring appears around the central peak representing the 
disordered structure observed in the inter-wall region, 
inside which the vortex density becomes more uniform 
as / approaches /=l/2 (compare (n v / c ) in (a) vs (b)). 
When the system is driven the ratchet induced periodic- 
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FIG. 4: Vortex structures: static b(n) obtained quenching 
from high T towards T—0 (L—128 but a fraction is shown for 
transparency. Vortices are white plaquettes). In (a) f—1/8 
and (b) /=3/4, both with (n v / c ) below, the number of vortices 
per column, (c) S(k)/f 2 corresponding to 6(n) in (a) with 
clear peaks due to the ratchet potential, as k\. (d) S(ki)/f 2 
vs I ac for /=l/8, 3/4 at T=0.05 and tj ac =10" 4 . 




FIG. 5: Time voltage response to an ac-drive of I ac =0.S for 
/=3/4 (smaller amplitude curve, black line), with (Vd c )>0, 
and 7 ac =0.7 (larger amplitude, red line), with (Vd c )<0. Insets 
(a),(b): instantaneous vorticity for 7 ac =0.3 and t/r ac — 0,0.25 
respectively, (c) Power spectrum for /=3/4 (solid black line), 
f—1/8 (dotted green line) and f—1/16 (dashed red line). 



ity survives but its amplitude decays with increasing I ac . 
We find that the degree of order decreases considerably 
for / around /=l/2 in the CR regime, although vor- 
tex walls at small I ac are still well defined, see FigfJJd). 
We now analyze the instantaneous vortex dynamics. In 
Figj5] we show examples of typical voltage response vs 
time during a period at which its steady value (Vdc) (see 
Figf2{a)) is already achieved. The amplitude of the al- 
most symmetrical temporal oscillations around Vt—0 is 
much larger, one order of magnitude, than (Vdc) for all 
I ac - This implies that rectification arises from tiny differ- 
ences between the two directions of motion and explains 
the low efficiency of the device and the need of very large 
T ac —ijJ~l to resolve the effect. The dc-limit analysis con- 
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firms this result [l9[- Insets (a),(b) in Figj5] show 6(n) 
snapshots at time [0,0.25]r ac respectively. Most of the 
time the moving vortex lattice is highly disordered (as 
seen at t/r ac — 0.25), except for small forces close to de- 
pinning (as in t/r ac ^0) where the ratchet-induced peri- 
odicity is more visible. This fluctuating order produces 
the small peaks in the averaged structure factor showed 
in FigfJJd). In Fig|5^c) voltage power spectrum for dif- 
ferent vortex densities are analyzed. Time correlations 
appear clearly at u ac and higher harmonics, as it is ex- 
pected for one particle in a ratchet potential. Note how- 
ever that interactions modify the harmonics structure, 
as it is seen comparing /=3/4 with more diluted vortex 
lattices, /=l/8 and /=1/16 where more harmonic peaks 
appear. We thus see that temporal order is only deter- 
mined by the ac-drive as it is expected for vortices moving 
incoherently. No trace of a washboard frequency is ob- 
served for any / nor I ac . Therefore no dynamical phase 
transitions are induced by the ac-drive and the system 
behaves as a modulated driven plasma. 

Discussion. Our results present some similarities 
but also some marked differences with other 2D vortex 
ratchet systems. First, let us note that inverted rectifi- 
cation appears only at low drives within a single range 
of high f around f=l/2 (FigHJ. This is similar to CR 
observed in arrays of dots in SC films by Villegas et al. 
0, but their explanation based on interstitial vortices 
moving over an inverted ratchet potential, is not suitable 
for JJA as interstitial can not be defined. Our results are 
also different to the multiple CRs observed and explained 
by de Souza Silva et al. in SC films with antidots arrays 
[7|. The differences can be attributed to the strong pin- 
ning generated by dots or antidots, which can localize a 
few vortices in very small regions and produce a very co- 
herent individual motion. In our case an extended weak 
potential pins and deforms collectively the vortex density 
although individual vortex motion remains highly inco- 
herent. In addition, underdamped dynamics was also 
shown to provide CR 0]. In our case no evidence of an 
effective dynamically induced inertia is found however, 
as no hysteresis in Fig[TJ nor out of phase voltage re- 
sponse is observed, FigfSJ CR in JJA thus presents novel 
features deserving discussion. 

Since CR is observed near the onset of dissipation, an 
analysis of the asymmetric depinning is useful. In the 
absence of vortex interactions it is clear that the criti- 
cal current in the easy direction of the ratchet potential, 
^dep V ~ AU/(p — 1) is smaller than the one in the hard 
direction ij^jj" ~ AU. Therefore the onset of rectifica- 
tion occurs when I ac ~ I^p V with a net motion in the 
easy direction. In presence of interactions the situation 
can be rather different since vortices tend to screen 
the ratchet potential thus modifying the critical currents 
and mobility of single vortex excitations in each direc- 
tion differently. It is clear that a small tilt of the ratchet 
potential in the hard direction produces well defined walls 



with a high vortex density at the local minima of the 
sawtooth. On the contrary a tilt in the easy direction 
favors a smoother vortex density profile. Calculating the 
force exerted on a vortex by a periodic array of vortex 
walls we can easily show that l\^ d is reduced by an 
amount proportional to the vortex density at the walls, 
Idep d= & U ~ & f wail, while in the easy direction I%£ y re- 
mains unchanged. Hence, CR can occur above a critical 
vortex density, f c , such that I^p d =I^ v at /o yield- 
ing fc ~ ^r(f5r)' This simple argument is consistent 
with the fact that inversion is observed at high /, dis- 
appears with increasing AU and is weakly dependent on 
p (details and dc limit confirming this scenario in [19(). 
At larger I ac the dynamical smoothening of the vortex 
density reduces the screening and the easy direction for 
rectification is recovered. Figfl] shows that this recov- 
ering occurs at higher I ac for higher /, consistent with 
the existence of denser walls at high /. Our results thus 
show that vortex-vortex interactions can overscreen the 
asymmetry of the ratchet potential as was noted in Q 
but the mechanism is purely collective, since individual 
vortex motion is highly incoherent like in a driven modu- 
lated fluid. Indeed, this effect could be relevant for a wide 
family of interacting many-body ratchet systems [l9j |. 
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